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The response of an unbounded and isothermal atmosphere to a point 
impulse of force is studied. Fronts of disturbances propagate with two 
different speeds corresponding to acoustic waves and to gravity waves. 

The distributions of kinetic energy density in the two wakes behind the 
two fronts are calculated. A striking result is that behind the acoustic 
front, the radiation density is very intense above and below the source, 
if the atmosphere is sufficiently stable to convection. This result 
may be of Interest in connection with the formation of spicule structures 
in the solar upper chromosphere* 
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1. INTRODUCTION 

In 4 previous paper (Kato 1965; hereafter this paper will be referred, 
to as I) we studied the response of an atmosphere to time-harmonic point 
disturbances. This atmosphere was assumed to be isothermal and unbounded 
and subjected to an uniform field of gravity. In the present paper we study, 
the response of the same kind of atmosphere to impulsive point disturbances, 
and emphasize some interesting effects of the degree of stability of the 

i 

atmosphere upon the response. 

The impulsive disturbance treated here will be assumed isotropic for 
two reasons, first, in paper I, we have already learned something "of the^_ 
directional response of the isothermal atmosphere by treating linear, 
oscillating forces aligned in different directions. In the case of impulsive 
disturbances we may expect that the differences in the direction will 
produce similar differences in the response of the atmosphere. Second, 
for an isotropic point impulse, the response of the atmosphere is spherically 
symmetric if the atmosphere is convectively neutral (7 * 1 ), as we shall see. 
Thus, we can concentrate our attention upon effects of the degree of 

stability of the atmosphere by studying deviation from spherical symmetry • 

« 

in the response. In the Appendix, however, we indicate how directed 
impulses may be treated and show how they may be combined to yield the 

results obtained with an isotropic source. / 

/ ; 

Now, Whitham (1961) showed that in the initial value problems of 

/ » 

propagation of waves , the amplitude of distances can be calculated easily 

■ , * •• • 

and directly with the help of some properties of group velocity. We apply 
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> ♦ , 

Whitham's method to study of the response of an atmosphere to an isotropic 

impulse, without solving an inhomogeneous wave equation. An advantage of 
such a direct argument is that it may bd taken over to a slightly non- 
isothermal atmosphere. 


II. GROUP VELOCITY 

' i 4 

In order to apply Whitham's approach, it is necessary to obtain an 
expression for the group velocity. If the dispersion relation is given by 
G(k. w) » 0, the group velocity, U, can be written as 




because 8 G/ 8 k^+ (9 g/8u) 9w/9k^- 0. In an isothermal atmosphere stratified 

under a' constant gravity, g ” (o, o, -g), the dispersion relation, G (k, w) » 0, 

is given by (e.g., see equation (5) in I) . • 

' . ... 

G(k, w) a - w 4 + u 2 (A 2 + « 2 ) - c 2 w 2 k 2 - 0 , (2) 

• < * 

■9 9 9 9 2 2 2 

where * k - k£ + k y + k^, k^ - k^ + k y , and 

• « x ■ r g/2 c, and w 2 - (7 - l ) 1 ^ 2 g/c , (3) 


• where 7 , g and c are the ratio of the specific heats, the magnitude of_ 
■ gravity, and the adiabatic sound velocity, respectively. 

If the angle between the direction of the group velocity and the 
positive z-axis is 9 and the longitude measured from x-axis is the' 
relation, U; H ~VG, shown in equation ( 1 ) can be written as 

2c 2 (w 2 - u 2 ) k 2 c 2 (w 2 -u 2 ) k 2c 2 w 2 k. 

z x z y z 


' v 2' x 
sin© cose 


sinG sin* 
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Expressing k for example, k and k by k with the use of equation (4), and 
Inserting them into equation (2) , we have an expression for k^ written by 
u, 9 ♦. Quite similarly we can have expressions for k^ and k^. These 
expressions axe l > 


/ 



u 4 (w 2 * -u 2 ) 


!(w 2 -u 2 )(w 2 -u 2 cos Z e). 


1/2 


slnB 



(5) 


and 


k “c 
z c 


(u 2 -« 2 )(u 2 -u 2 ) 

2 2 2 
(« -cCJeo e 9) 


1/2 


CO80 


( 6 ) 


Inserting' the values for k , k and k given in equations (5) and (6) 

x y z 

into a component of equation (1) , the magnitude of the group velocity can 

i 

be expressed as (see also Moore and Spiegel 1964) 


- c 2 


> 


(u 2 -w 2 )-\u 2 -uj) (u 2 -w 2 cos 2 e) 
u 2 [(u 2 -u 2 ) 2 + u 2 (u 2 -w 2 )sin 2 e] 2 


(7) 


Of course, the group, velocity has meaning only in the frequency ranges of 

•a 

acoustic waves (» > «0 and of gravity waves (« < » 2 cos 9). Here we 

should mention that in waves radiating energy in the critical direction, 

0 (cos 2 0 - u 2 /u 2 ) , of gravity waves the group velocity is zero (see 
c c Z 

equation (7)) and the wavelength is infinite (see equations (5) and (6)). 


III. FREQUENCY DISTRIBUTION IN WAKES 

We assume a point Impulse of force at the origin, * r ** 0» at the 
time, t - 0. As the medium is dispersive, there remain wakes after the 


1) We consider the plane wave exp(i«t-ik.r) for « > 0. If we consider the 

plane wave exp(iwt+ik.r) for u > 0 the signs of k^, k^, and k^ in 

equations (5) and (6) should be changed. 


rrjo 
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'passage of fronts of disturbances. At points in the wakes far from the 
fronts the disturbances are broken down into wave trains, the different 
wave number components or frequency components distribute over the physical 
jr-space, and each k (or u) propagates with the group velocity (Whithaxn 1961 

see also Appendix B in Lighthill's paper 1960). That is, far from the 

• • 

fronts of disturbances , k (or «) and x are related through the group 
velocity U by 


r <■ U (u, 6) t • (8) 


This relation enables us to calculate the frequency distribution in wakes 

except for the region near the fronts. We introduce a sphere of radius 

r ** ct (this sphere represents the front of aeons tiq waves). Then, from 

equation (8), we have r/r * U(w,0)/c or by solving this for w we obtain 

c 

u a u(r/r .0:7) , which shows the frequency distribution in wakes. In 
• c 

t' ' • 

general, u is not a single valued function because there are two wave modes. 

• i ■ * . 

(the acoustic and the gravity wave modes). The surfaces of constant values 
of w in the ^r-space are drawn in Figures 1 and 2 (see also Figure 4 in 

r ' • 

Moore and Spiegel's paper (1964)) for the acoustic wave mode (« > «^) 

and for the gravity wave mode (« < « 2 )» respectively. The three curves 

for each value of « in Figure 1 are for 7 ■ 1, 4/3 and 5/3 respectively. 

The left hand side in Figure 2 is for 7 ■ 4/3, and the right handi side is 

for 7 " 5/3. • 

Figure 1 shows that a front of acoustic waves propagates spherically. 

1/2 

from the source with the sound velocity (7 R T^) . After the passage of 

the front there remains a wake. In the wake frequencies of disturbances 
tend rapidly to the upper characteristic frequency of the atmosphere, w i»' 
as shown in Figure 1. The result that there is a wake behind the acoustic 


7 


front, and the frequencies of disturbances in the wake tend to the upper 
characteristic frequency is qualitatively the same as in the one-dimensional 
case studied by T-anih (1908). (Lamb studied the vertical propagation of 
disturbances from an impulsive force acting uniformly on an infinite 
horizontal plane.) The other front (front of gravity waves) is the one 
formed by the envelope of disturbances propagating from the source as . 

gravity waves, and after the passage of the front, wave motions at a 
point in the wake consist of two different frequencies as shown in Figure 2. 
The gravity waves, of' course, do not appear in Lamb's one-dimensional case. ' 
Here we note that Figures 1 and 2 can also be interpreted as polar 
diagrams of magnitude of the group velocity for acoustic waves and for 

i 

gravity waves, respectively. That is, the radial distance from the origin 
of a point on a curve in Figures 1 and 2 denotes the magnitude of the group, 
velocity in. that direction and for that frequency. - — . 

IV. DISTRIBUTION OF KINETIC ENERGY 
DENSITY IN WAKES 

a) Whitham's Method 

Whitham (1961) has shown that in the initial value problem the amplitude 

variation of disturbances can be calculated easily and directly with the 

help of some properties of group velocity. We shall apply here Whitham's 
» 

method to' calculations of the distribution of kinetic energy density in 
wakes radiated from an isotropic impulse. Following Whitham, we consider 
plane wave solutions, exp(i«t - ik.r) , where « and k are related by 
the dispersion relation, G(k, «) - 0. In our present problem, the dis- 
persion relation is given by equation (2) . This dispersion relation shows 


that for a given Jk there are always two) positive values of « Z ; one corresponds 
to the acoustic wave mode and the other to the gravity wave mode. For 
disturbances associated with the each wave mode, the magnitude of velocity, 

j. 

u, times' one-half of the undisturbed density, P Q » may be analyzed into 
plane waves as 


P * u(r 
o ~ 


oo 

(r,t) -J F(k) exp(iut - ik-r) dk, (w^ 0) 


where « is related to k^ by the dispersion relation, where F can be ex- 

t . * 

pressed in terms of the initial conditions by 


>!v> I 


Pq »(r,o) exp(ik-r) dr, . 


The right hand side of equation (9) is complex and its real part gives the 
magnitude of the disturbance. 

If an impulse is isotropic' the form of the initial condition is 
obvious, and we can apply Whitham' s method easily. That is, the initial 
condition can be written as 


• P Q u(r,o) - A 6(r) , A - complex constant. 

Thus, inserting equation (11) into equation (10), we have 


*<& — ^3 
<2 «) J 


That is, F(k) depends neither on the direction of k nor on .its magnitude. 
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The total kinetic energy of disturbances at the time, t, is given by 


J 1 p o % u( 5' c) p 0 % u * ( £’ c) d £ * 


(13) 


and, on substitution of equation (9), this can also be written as 

» OP 

. (2«) 3 J \ F(k) F*(k) dk . (14) 

As mentioned in the previous section, after the initial disturbances have 
been broken down into wave trains the different values of k are distributed 
over the physical r-space and each value of k propagates in the physical 

i 

r-space with the group velocity. Therefore, we have a direct correspondence 
between equations (13) and (14) . That is , if the wave number at position 
r is k and the volume element Ar contains waves with wave numbers in an 
element of wave number space, then the contribution of that element to 
energy is . 


\ P q ^ u(r, t) u* (r, t).A r , 


(15) 


or in wave number space. 


| (2m) 3 F(k) F* (k) Ak 


(16) 


and the two must be equal. Therefore the amplitude of p • u (r, t) is 




(17) 


If we follow the wave number k and retain a constant Ak, the volume ele m e n t 
Ac will vary due to the small differences in the group velocity for the 
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different wave numbers in the element. Using this fact, we can calculate. 

Inserting the expression for Ak/£r calculated in this way into 
equation (17), Whitham (1961) obtained 


“<£>‘>1 * ( t ) ' l F ®l 


(18) 


The above is a direct and simple derivation of the intensity of disturbances 
(compared with the method of stationary phase approximation (Llghthill, I960)) 
although the derived information is restricted. (Equation' (18) shows only 
the magnitude of velocity, and components or the phase of the velocity are 
not obtained.) 


b) Application to an Isothermal Atmosphere 

To apply , Whitham' s method to an atmosphere we must calculate 
2 

det 9 <*>/9kj8kj explicitly as a function of u. Considering that in an 
isothermal atmosphere the dispersion relation, G(k, u) ■ o, is given by 
equation (2) and also using equations (5) to (7), after rather lengthy 
calculations we have the following relation 


det 


8 2 <o 




c: 

8 


(m 2 - 


2 2 2 
cos 0) 


w 3 (« 2 -u 2 ) 


U 2 ^ , («£<>). (19) 
du 


With the use of equation (19), the amplitude of p * u(r, t) given 

o ^ 

by equation (18) can be written as 


«(£,t)|- 


2jc 

t 


3/2 , 

l*Q0| 


« 3 (w 2 - CO 2 ) l 

L (co 2 -u 2 cos 2 e ) 2 U 2 8U 2 /8« 2 J 


( 20 ) 


For an isotropic impulse F^k) is constant and given by equation (12). 
Consequently, the kinetic energy density Eg Or, t), associated with a 


wave mode is 
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\ i | »o % 2 


(2it) 3 C 2 t 3 


12 a 3 (J- 4) i_ (B>0) 

(u-u^cose) iTdir/Bur / 


(cf. equation (36) in the Appendix) where w is the frequency of the waves 

• • t 

at the position r and the time t and is given, by equation (8) as discussed 
in Section III. 

t t 

Here it should be remembered that there are three cases where equation (8) 
has zero, one and two real positive solutions of « for a given jc and. t (see 
Figures 1 and 2). These three cases respectively correspond to the three 
cases; i) where, the front of acoustic waves has not reached the observation 
point (r , t) yet; ii) where the front of acoustic waves has passed the 
'observation point, but the front of gravity waves has not reached there yet, 
and Hi) where both fronts have passed. Equation (21) is valid for each 
wave mode. Thus, in case ill), the total kinetic energy density, at jr and t 
is given by summing equation (21) for the acoustic wave and equation (21) 
for. the gravity wave (cf. equation (365 in the Appendix). 

Finally, it should be mentioned that the attenuation of intensity like 

—3 —n « 

t or r (r id the n-dimensional space) along any radius vector is 

characteristic of impulses propagated three dimensionally outward in a 

dispersive medium. An impulse contains a range of frequencies ,. so 

3 

’ disturbances propagate as a volume and their volume increases as r . 


V. NUMERICAL RESULTS FOR INTENSITY DISTRIBUTION IN WAKES 


Equation (21) can be rewritten as 


\ “ 3~3 

* ~ (2j 0 V 


U ) 3 

Ul 2 .^- 


F («,e;7), (22) 



(23) 

and C ± denotes the iso thermal sound velocity. The coefficient of F in 
equation (22) is independent of both the ratio of specific heats, 7, and 
position in the ^r-space. The space and 7 dependencies of (r, t) are 
represented by the non-dimensional quantity, F(w,0;7), defined by 
equation (23). We will calculate F(m, 0;7) due to acoustic waves and 
gravity waves, separately. 

% •. * 

a) Behind the Front of Acoustic Waves 
As discussed in Section III, the distribution of <•> in the r-space 
is a function of r/r , 0 and 7. So, F is also a function of Tr/r , 0 and 7. 
The calculated distribution of F(«, 0; 7) in the jr-space is shown in 
Figures 3(a), 3(b) find 3(c) for cases of 7 “ 1, 4/3 and 5/3, respectively. • 
Here we mention that we do not calculate the intensity of disturbances 
near the front because near the front disturbances do not break down into 
wave trains yet and the approximation used to obtain equation (21) is not 
applicable. That is, a point in the k-space does not correspond to a 
point in the ^r-space exactly and energy within a given volume in the k-space 
is not conserved. 

1 

Figure 3(a) shows that intensity of disturbances in the wake decreases 
monotonically and spherically after the passage of the front in a con- 
vectively neutral atmosphere (7=1). This is natural because~an isotropic 
impulse acts in an atmosphere having no particular direction (convective 
neutrality).. However, figures 3(b) and 3(c) show that after the passage 
of the> front radiation density is very intense above and below the source 


■ where 


3 2 2 

w J (oT-tO 

F(w,0*7) “ 3 2 2 2 2 

uj(u-ttjcos0) 


— % 4 

7 * C* 

u^i^/au 2 


Cor 
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if .the degree of stability of the atmosphere increases (i.e. if the 

value of 7 increases). This result is suggested by the following two 

fatts. First,, an impulse contains all frequency components. Second, 

* 

the results of paper I show that the atmosphere can respond strongly to 
a time-harmonic disturbance whose frequency is near the upper characteristic 
frequency, u^, in a sharply restricted region above and below the source, 
if the atmosphere is sufficiently stable to convection (see Figure 1(a) and 
Figure 3 in paper I). 

This striking result shown in Figure 3(b) and 3(c) (radiation density 
is very intense sharply above the source) may be of interest here, because . 

, i 

this may suggest a reason of formation of spicule structures in the solar 
upper chromosphere. That is, vertical focussing is a. possible result of 
stratification alone. 



b) Behind the Front of Gravity Waves 

The distribution of F(w, 9; 7 ) behind the front of gravity waves are 

* calculated for the case of 7 * 5/3, and added in Figure 2. On the front,. 

9U /3w is zero because the front is defined as the envelope of disturbances 

of all frequency components. So, if equation (23) is applied literally, 

intensity becomes infinite on the front. Equation (23) is, however, invalid 

near the front for the same reason mentioned in the previous sub-section 

for the case of the front of acoustic waves. 

Here we should mention that intensity shown in Figure 2 comes mainly 

2 2 2 

from waves radiated near the critical direction 0 c (cos 9 C * “ / u 2 ^* 

v 

As mentioned before (Section II or paper I), the wave number of waves 
radiating energy in the critical direction is infinite. Thus, a careful 
discussion whether intensities of gravity waves shown in Figure 2 can be 



OOP - 


1 . 4 . 

radiated to a great distance from the source in actual physical situations 

' » • 

may be necessary. 

In conclusion, the author would like to express his sincere thanks to 

I* 
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Institution, 1964. The author would like to thank Prof. Willem V. R. Malkus 
and other participants for discussions. He would like to thank the National 
Academy of Sciences for its support and Dr. R. Jastrow for extending to him 

i 

the hospitality of the Goddard Institute for Space Studies. 
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* APPENDIX: RESPONSE TO A POINT IMPULSE ACTING IN A PARTICULAR DIRECTION 

» ' 

a) The Asymptotic Solution of the Inhomogeneous Wave Equation 
As. done in this paper, the response of the atmosphere to an impulse 
can be studied as an initial value problem. However, the same problem 
can be also studied by solving the inhomogeneous wave equation. The 
. later method vill be simpler in the case where the impulse is a directed 
point source. Let us asstime an external force, f (x, t), in an unbounded 

i , 

and isothermal atmosphere stratified under a constant gravity. Then, the 
inhomogeneous wave equation represented by the pressure variation, p^, 
is (see equations (5) and (6) in 1) 


+ ~Z <° 2 V 2 “ -?> + * «2 V? 

8 c B t Z L ZVi 


(p l /p o } “ S » (24) 


where 


t) - - c 2 ^ , 8z 


[h + “2 /s 


f z <£’ 


% 2 
- P D C 


3* 


+ u 2 — f (x, t) + — f ( X, t) 
L 0t“ 2 J19 x x ~ 8y 7 V 


(25) 


We solve equation (24) under the condition that dependent variable, 
p^/p o , is everywhere zero until an initial instant, t =* o, and at the time, 
t = o the force, f, acts instantly. 

We will introduce a four dimensional space (r, ct), and a four 
dimensional wave vector space (k, «/ c). Then, using the approximation of 
stationary phase as t and r' *•*■* (whefe r' is the distance between the 
observation point, r, and the source region x ; hereafter we take the 
origin in the source region and do not distinguish r' and r because r' *♦ °°) , 
the asymptotic solution of equation (24) is (Lighthill, 1960) 



% ^ ff , 12*L 5/2 C o *<*» w/c > «pU(-Jc-i+ «t)l 


VVfe*» fI "^372 


QgI • Ik I 1/2 
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» (26) 


^ 2 2 2 1/2 
where r is the four dimensional distance, (r + c t ) , and where the 

sum 2 Is over all points (u^o) on the hyper-surface, G(k, «/c) - 0 (see 
equation (2)), whose normal Is anti-parallel to the four dimensional obser- 
vation point (r, ct); that Is, over all points where 


r : ct -VGOc, «/c) ; - G (k , «/c) . (27) 

(cf. equations (1) and (8)). Moreover, in equation (26), u/c) is a 
Fourier component of S(x, t); 

®(lc, w/c) - — 1 -t f S(x, t) exp[-i(-k-x + «t)] dx d(ct), (28) 

~ (2 ity J ~ 

C Q is a phase factor of modulus 1; OG Is the gradient of G in Cartesian 
four space; and K ,is the Gaussian curvature of the hypersurface - each 
element of the matrix of products of first derivatives of G being 
multiplied into the corresponding co-factor of the matrix of second 
derivatives, and divided by G 


b) Response in Kinetic Energy 


As a measure of the response we adopt kinetic energy dentlsy, E^Cr, t) 
Kinetic energy density is related to pressure fluctuation by an equation 
(see equation (18) in paper I). Inserting equation (26) into the relation, 
after a lengthy calculation we have 




^(r, t) 


12*1 

4t 3 


(«^) J 


C 7 (u 2 *u 2 ) 3/2 (u 2 -w 2 cos 2 6) 1/2 


V l®l' 


U9U 2 /8m 2 


(29) 


vhere the prime added to Z means that the summation is performed only for 
points vhere u >-o of the points satisfying equation (28), 

In deriving equation (29) we used 

.|OG| 2 -4C 2 « 4 ■ < - p~^ > (l + 4) (30) 

(u *-uqcoa B) ‘ 1 r • 


and 




Idol 


2 2 2 3 2 2 5 

u (u -up J (w -U^) 

L (u 2 -u 2 cos 2 0) 


1/2 


1 8P 2 
U 3 8w 2 


(31) 


which are obtained from the definition of G end K and the relation be* 

, * »' , , 

tween (k, w/c) and. (r, ct) shown in equation (27), 

As 


' c) Response to Three Kinds of Impulses 

To calculate I® I 2 in equation (29) explicitly, we must specify the 
external force, f. We study about the following three kinds of sources, 
separately; that is, point Impulses acting to x, y and z directions at. 


the origin, r ■ 0, and at 

the time. 

t - 0: 

i 


• • •’ 1 

M 

| 

[M.' 

case (a) 


f i“ : 1 

i 

6(x)5{t) 

AS ’ 

M 

case (b) 

(32) 

•: v 1 



UJ 

case (c) 

« 

• 
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In case (a) , for example, Inserting equation (32) into equation (25) and 

I 

.i * 

performing Integrations with respect to x and t to obtain 9 defined by 
equation (28), then we have 0 expressed by k and fa). We eliminate k from 
this egression for 9 with the help of equation (27) and square the 
expression. Then, we have 



2 2 2 2 2 

*1 -14 2 2 

h c “ — 5 T *“5 5 sin 0 cos ♦, case (a) / 

(2*) 8 00 (fa> 2 -fa)^cos 2 e) 

( 33 ) 


•where the suffix 1 of 0 denotes case (a), and where p denotes the 

oo 

undisturbed density at the position of the impulse. Quite similarly, for 

I . 

cases (b) and (c). We have 


a 2 (fa) 2 - « 2 )(fa) 2 -fa> 2 ) 

1 9-1 ^ u* — , — s = sin 2 0 sin 2 * , case (b) 

2 (2*) 8 00 (fa) 2 -fa) 2 cos 2 0 


(34) 


and 


l9 3 l ! 


4 


-1 > .4 


( 2 *) 


8 P oo 


2 .2, 


C- faT [(upu)p + 


(fa) 2 -fa) 2 ) (u 2 -w 2 ) 

2 2 2 
(fa) -u^cos S) 


2 

cos 0 ] 


case (c) 
(35) 


If equations (33), (34) or (35) are inserted into equation (29), we have 
expressions for distribution of kinetic energy in the ^r-space resulting 
from impulses directed in the x, y, or z directions, respectively. 


d) A Relation to the Responce to an Isotropic Impulse 
To show a relation between the response of the atmosphere to an 
isotropic Impulse and to the directed impulses , we assume three atmospheres 
of the same kind, and assume that in each atmosphere an impulse acts in 
the x, y and z-dlrection with a same amplitude (A q 5 * A^) at the 
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same point and time, respectively. We sum the kinetic energy density at 
the corresponding points in these three systems at the same time. That is, 
we sum equations (33) to (35) and insert the sum into equation (29) to 
obtain 


Ve.‘) 


i 

3 3 
4(2rt)V 





(w 2 -u 2 cos 2 ©) 2 U 2 3U 2 /9 u 2 



As anticipated, the intensity distribution of the kinetic energy density 
in the ^r-space shown in equation (36) is quite same with that for an 
isotropic point impulse shown in equation (21). 
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Fig. 1 


Fig. 2 


Fig. 


« < 

Figure Captions 

t 

. Frequency distribution in wave motions in the wake of acoustic 

waves, r denotes the position of the front of disturbances, 
c 

Three cases (7 ” 1, 4/3 and 5/3) are drawn on the same" map. 

. Frequency distribution in wave motions in the wake of gravity 
waves. The envelopes denote the front of disturbances. The 
left hand side is for y ■= 4/3; the right hand side is for 
7 ** 5/3. The broken lines denote intensity of disturbances 
radiated by an isotropic point impulse, and the numerical values 
represent values of F(w, 6 ; 7 ). The circular arcs at both sides 

1 

of the figure denote the position of the front of acoustic waves. 

» Distribution of kinetic energy density in the wake of acoustic 

* / 

waves radiated by an isotropic point impulse. The numerical/ 

values added to iso-intensity curves represent values of 

F(w, ©; 7 ) defined by equation (23). r denotes the position 

c 

of the front of disturbances. Figures 3(a), 3(b) and 3(c) 
correspond. to the cases 7 ■ 1, 4/3 and 5/3, respectively. 


